arXiv: 1504.05545v2 [hep-lat] 3 Nov 2015 


Lattice QCD Green’s functions in maximally Abelian gauge: infrared Abelian 

dominance and the quark sector 

Mario SchrodO 

Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Roma Tre, 00146 Rome, Italy 

Hannes Vogl[] 

Institut fur Theoretische Physik, Auf der Morqenstelle 14, 72076 Tiibinqen, Germany 

(Dated: November 4, 2015) 

On lattice gauge field configurations with 2 + 1 dynamical quark flavors, we investigate the mo¬ 
mentum space quark and gluon propagators in the combined maximally Abelian plus 17(1)3 XU( 1)8 
Landau gauge. We extract the gluon fields from the lattice link variables and study the diagonal and 
off-diagonal gluon propagators. We find that the infrared region of the transverse diagonal gluon 
propagator is strongly enhanced compared to the off-diagonal propagator. The Dirac operator from 
the Asqtad action is inverted on the diagonal and off-diagonal gluon backgrounds separately. In 
agreement with the hypothesis of infrared Abelian dominance, we find that the off-diagonal gluon 
background hardly gives rise to any non-trivial quark dynamics while the quark propagator from 
the diagonal gluon background closely resembles its Landau gauge counterpart. 
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I. INTRODUCTION 

Despite the long standing acceptance of quantum chro- 
modynamics (QCD) as the correct theory to describe 
the strong interactions of quarks and gluons, a basic un¬ 
derstanding of its main characteristic features, the dy¬ 
namical breaking of the chiral symmetry and confine¬ 
ment, is still lacking. Several scenarios for the underlying 
mechanism of confinement have been suggested over the 
decades. Under the most popular ones are the Kugo- 
Ojima [li and the Gribov-Zwanziger scenarios HE] , as 
well as the dual superconductor picture HI]. Some com¬ 
mon aspects of the different confinement criteria have 
been investigated in um- 

The dual superconductor picture of confinement is es¬ 
pecially appealing since it offers a rather intuitive ap¬ 
proach to confinement. Type II superconductors in their 
superconducting phase are known to repel external mag¬ 
netic fields below a critical value of the external field 
strength. If the magnetic field exceeds that value, tubes 
of magnetic flux (Abrikosov vortices) begin to penetrate 
the superconductor. The flux tubes are encircled by 
Cooper pairs which squeeze the latter. Identifying the 
penetrating magnetic flux tubes with the color electric 
field of the Yang- Mills vacuum, and moreover the con¬ 
densed electric monopoles (Cooper pairs) with color mag¬ 
netic monopoles, one finds a dual picture where hypothet¬ 
ical magnetic monopoles at the beginning and end of the 
Abrikosov vortices correspond to the confined quarks of 
QCD. Therefore, it is suggested that confinement is due 
to the condensation of color magnetic monopoles which 
squeeze the color electric flux tube between quarks and 


* mario.schroeck@roma3.infn.it 
t hannes.vogt@uni-tuebingen.de 


antiquarks. The dual superconductor picture has a far 
reaching consequence: the Abelian parts of the gauge 
fields should dominate the nonperturbative infrared (IR) 
dynamics [12] . 

While confinement is the reason why individual quarks 
and gluons have never been observed in experiment, the 
theory still allows to investigate correlation functions of 
single quark and gluon entities: the QCD Green’s func¬ 
tions. QCD is a gauge theory and therefore the gauge has 
to be fixed in order to study the fundamental two-point 
functions. The maximally Abelian gauge, as the name 
suggests, rotates the gauge fields such that the diago¬ 
nal, i.e., Abelian parts of the gauge fields are enhanced 
over the off-diagonal parts. This renders the maximally 
Abelian gauge particularly suitable to study IR Abelian 
dominance. 

Lattice QCD provides an approximation to the con¬ 
tinuum formulation with a finite number of degrees of 
freedom which allows to perform numerical simulations. 
Various attempts of investigating the dual superconduc¬ 
tor picture in lattice QCD have been carried out, see, e.g., 
muni. IR Abelian dominance has been demonstrated 
in lattice gauge field theory in SU(2) fT5HT7] and more 
recently in SU(3) [TS1 ITT?] by studying the infrared be¬ 
havior of the gluon propagator. In the recent study m 
almost perfect Abelian dominance of the string tension 
on large physical volumes in quenched SU(3) has been 
found. Complementary to the lattice approach, Abelian 
dominance has been analyzed perturbatively and 

in Dyson-Schwinger and renormalization group equation 
studies [251125] . 

It is interesting to study the influence of Abelian domi¬ 
nance on chiral symmetry breaking and the hadron spec¬ 
trum; in [SB] it has been found that SU(2) Abelian pro¬ 
jected fields give a chiral condensate which closely re¬ 
sembles the results of strongly coupled gauge theory. In 
m quenched hadron spectra in Abelian gauge fields, ex- 
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tracted by maximal Abelian projection have been stud¬ 
ied: the ratios of the hadron mass to the square root 
of the string tension of the Abelian fields are similar to 
those of the full SU(3) theory. The authors concluded 
that they have found Abelian dominance (and monopole 
dominance) for the hadron spectra. 

In the current paper we advance previous investiga¬ 
tions of the gluon propagator in maximally Abelian gauge 
from pure Yang-Mills theory nil [191 to full dynamical 
QCD by adopting Nj = 2 + 1 gauge field configurations. 
Furthermore, we analyze for the first time the maximally 
Abelian gauge quark propagator. In order to obtain in¬ 
sights in the dependence of chiral symmetry breaking and 
the dynamic mass generation of quarks on the type of 
gluon background, we invert the Dirac matrix separately 
on the diagonal and off-diagonal gluon fields. 

The remainder of this work is structured as follows. 
In Sec.jTTJwe introduce the maximally Abelian gauge and 
discuss some aspects of its implementation on the lattice. 
After reviewing the methods to extract the QCD Green’s 
functions on the lattice in Sec. m we list details of our 
lattice setup in Sec. |IV|and present our results. 

II. THE MAXIMALLY ABELIAN GAUGE ON 
THE LATTICE 

The continuum gauge fields are given by 

A^x) = l -j2\Af{x), ( 1 ) 

Z i= 1 

where the \ are the Gell-Mann matrices and ( x ) are 

real. On the lattice, the latter translate to the lattice 
link variables U^(x) £ SU(3). The continuum and lattice 
fields are related via 

U„ (x) = e iaaoA ^ x) (2) 

with a being the lattice spacing and go the bare coupling 
constant. A gauge transformation in the language of lat¬ 
tice QCD reads 

Up{x) -* g(x) U^(x)g(x + A) 1 (3) 

with local gauge transformations g(x) £ SU(3). 

The maximally Abelian gauge (MAG) aims at min¬ 
imizing the off-diagonal part of the gauge fields, i.e., 
(x) with i ^ 3,8. This is equivalent to maximizing 
the following functional of the link variables, 

^MAG3 [^1 = E tr [W*) WX)1] (4) 

X,fl 

+ tr [X s U l _ l (x)X 8 U IJ ,(x^] (5) 

where A 3 and X 8 build the Cartan subalgebra of SU(3). 
Once the functional Eq. Q resides in a local maximum, 


the gauge condition 

y^E(E u ^ ) ( x )°' 3U M ) ( a: ) t ( 6 ) 

d x,j ll 

+u ( j?\x - A) f (Tzvffix - £)) (7) 

becomes small. Here N,i is the number of Euclidean 
spacetime indices, V the number of lattice sites and the 
u^\x), j = 1,2,3 are the SU(2) subgroup elements of 
the link variables £7/x(x)rl In practice we reach a gauge 
precision of 9 < 10~ 13 . More details of the implementa¬ 
tion can be found in Refs. [28’, 29]. 



FIG. 1. The relative frequency of the gluon field components 
A^\x) (in lattice units) from a single gauge field configura¬ 
tion in Landau gauge and MAG. In Landau gauge the eight 
field components follow the same normal distribution while 
in MAG a trend of the off-diagonal components (MAG-Ojff) 
towards smaller values compared to the diagonal components 
(MAG-diag) is manifest. 

Once a maximum of Eq. 0 has been reached, the di¬ 
agonal gluon fields, A^\x) with i = 3,8 are favoured 
over the off-diagonal gluons. In Fig. [T] we show the dis¬ 
tribution of the gluon fields A^}{x) from a single lattice 
gauge field configuration. While in Landau gauge none 
of the gauge fields is favoured and thus their distribution 
lies on top of each other, the distribution of the MAG 
fields shows a clear shift of the diagonal fields towards 
larger values and correspondingly a shift towards smaller 
values for the off-diagonal parts. 

The MAG functional Eq. Q is invariant under gauge 
transformations of the form 

gd{x ) = exp A 3 -Mo/^Asj (8) 

and therefore the MAG is an incomplete gauge condi¬ 
tion; it leaves a remaining U( 1) 3 x XJ{ l)s gauge freedom. 


1 Two SU(2) matrices overlap on the diagonal of U^(x) and the 

third one consists of the corners of 
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The latter we remove by enforcing that U^x) fulfils in 
addition to the MAG the Landau gauge condition 

^La„dauM=^E tr \g(x) U^{x) g{x + A) T ] —■> max., 

/X ,X 

(9) 

with respect to diagonal gauge transformations g(x) = 
gd(x). 

Lastly, in order to study the diagonal and off-diagonal 
parts of the gluon fields separately, we extract the fields 
A^x) from the link variables Eq. ([2]) via the first order 
approximation 

A^x) = (U^x) - C/ M (x) t ) . (10) 

Ziag o traceless 

Note that we do not make use of the “exact” logarith¬ 
mic definition of the lattice gluon fields in order to stay 
consistent with the definition of the gluon fields in the 
Landau gauge condition [28] {^] 


For the MAG case we split the propagator in a diagonal 
^l ag (fc 2 ) = \ E (A«(fc)AW(-fc)) (15) 

z=3,8 

and an off-diagonal part 

D f^ k2 ) = l E (4 l) (fc)4°(-fc)) ■ (is) 

i^3,8 

Due to the residual 17(1)3 x t/(l) 8 Landau gauge fixing, 
the diagonal propagator is transverse, whereas the off- 
diagonal propagator has a longitudinal and a transverse 
component 



B. Quark propagator 


III. QCD PROPAGATORS 
A. Gluon propagator 

In Landau gauge the gluon propagator in momentum 
space is transverse and diagonal in color space 

D^(k 2 ) = (A®(k) A P(-k)) ( 11 ) 

° 2=1 

= (V-^)^ 2 ), ( 12 ) 

where A^\k) are the Fourier transformed gauge fields 
extracted from the links by Eq. ( fl0| ). Appropriate for the 
Symanzik-improved Liischer-Weisz gauge action [3T], we 
define the momentum variable as 



are the discrete lattice momenta. The transversality of 
the momentum space propagator is a direct consequence 
of the Landau gauge condition d^A^x) = 0. 


2 We performed some checks using the logarithmic definition: the 
qualitative behavior of the gluon propagator is the same, as it was 
also found in Ref. m■ The main difference here is, as expected, 
that the diagonal part of the longitudinal gluon propagator from 
the logarithmic definition does not vanish. 


In manifestly covariant gauges, the interacting quark 
propagator S(p\p 2 ), renormalized at the renormalization 
point /i, can be decomposed into Dirac scalar and vector 
parts 

S{p-,p 2 ) = (i]/)A(ij,;p 2 ) + B(fj,-p 2 )) 1 (18) 

or equivalently as 

S(p',p 2 ) = Z(p-p 2 ) {ij/> + M(p 2 )) 1 . (19) 

In the last equation the wave-function renormalization 
function Z(p\p 2 ) = 1 /A(p\p 2 ) carries all the information 
about the renormalization scale and the mass function 
M(p 2 ) = B(g;p 2 )/A(p-,p 2 ) is a renormalization group 
invariant. 

The lattice regularized quark propagator Sl{p 2 ',ci ), 
which depends on the lattice spacing a, can then be renor¬ 
malized at renormalization scale /i with the momentum 
independent quark wave-function renormalization con¬ 
stant Z 2 (g;a), 

S L {p 2 ;a) = Z 2 (p,;a)S(p-,p 2 ). (20) 

The momentum subtraction scheme (MOM) has the 
renormalization point boundary conditions Z(p; p 2 ) = 1 
and M(p 2 ) = m(p) where m(g) is the running mass. 

The nonperturbative functions M(p 2 ) and Z{p 2 ) = 
Z 2 (p\a)Z{p]p 2 ) can be extracted directly from the lat¬ 
tice. To this end we invert the Asqtad fermion matrix (32) 
in order to obtain the quark propagator which we sub¬ 
sequently Fourier transform to momentum space. Tak¬ 
ing basic Clifford algebra properties into account we can 
extract the dressing functions. For details we refer to 
(33J[3J. Note that the lattice dressing functions will be 
functions of the lattice quark momenta (which differ from 
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the gluon momenta (13)) and for the Asqtad action these 
are defined by 


K = sin(^) ( 1 + - sin 2 (p M ) ) . (21) 


We perform a cylinder-cut [35] on all our data and aver¬ 
age over the discrete rotational and parity symmetries of 
Sl(p 2 ',ci) to increase statistics. 


IV. RESULTS 

A. Gauge configurations 

For our simulation we adopted two sets of gauge field 
configurations generated by the MILC collaboration [ 13 - 
[55] : a “coarse” set of size 20 3 x 64 with lattice spac¬ 
ing a = 0.12 fm, which consists of five dynamical plus 
a quenched ensemble, and furthermore a “fine” set con¬ 
sisting of a single ensemble of size 40 3 x 96 with lattice 
spacing a = 0.09 fm. The configurations were generated 
with the Symanzik-improved Liischer-Weisz gauge action 
[40] and have been made available to the lattice commu¬ 
nity via the Gauge Connection BH Both sets include 
two light degenerate (Z) and one heavier quark flavor (s) 
(except the quenched ensemble), implemented with the 
Asqtad improved action [35]. The parameters of the lat¬ 
tices are summarized in Table [T] for the reported lattice 
scales and quark masses we refer to the original work 
[36H39]. 



TABLE I. Overview of the gaugefield parameters: the lattice 
size IV 3 x Nt , lattice spacing a, dynamical quark masses mi 
and m a (oo indicating quenched gauge fields) and the number 
of configurations that enter our analysis. 


propagators (diagonal transverse, off-diagonal longitudi¬ 
nal and off-diagonal transverse) we include the Landau 
gauge propagator for comparison. The Landau gauge 
propagator has been studied on the same dataset in [35] . 
Following their setup we use the same renormalization 
condition 


D(k 2 = p 2 ) = \ (22) 

at p = 4GeV. In the IR we find suppression of the 
propagator with dynamical quarks due to screening ef¬ 
fects. Compared to Landau gauge, this effect is more 
pronounced in the maximally Abelian gauge. Decreasing 
the quark mass leads to a further suppression in the IR, 
however the dependence on the quark mass is small. 

In Fig. [3] we present the gluon form factors of the fine 
ensemble. It is very important to note that in both 
Figs. [2] and [3] it is obvious that the diagonal parts of 
the MAG gluon propagator are pronounced as compared 
to the Landau gauge counterparts, and, respectively, the 
off-diagonal parts are suppressed. This extends previous 
findings of IR Abelian dominance from pure Yang-Mills 
theory to full QCD. 

Our findings for the fine ensemble in Fig. [3] are in very 
good agreement with the SU(2) results of [H] considering 
the fact that we study the SU(3) propagator and included 
dynamical quarks. There the authors found a wide max¬ 
imum of the dressing function k 2 D^f(k 2 ) at around 2 
GeV and a sharp peak of k 2 D dlag (k 2 ) around 0.7 GeV. 
To compare our results with the quenched SU(3) results 
of [T 93 we applied a fit to our data with the their function 


D(k 2 ) 


Z 

{k 2 + m 2 )" 


(23) 


in the same momentum regime k < 3 GeV. Qualitatively, 
the results compare well to the quenched results from 
Ref. [T5] . 



m [GeV] 

V 

Z 

X 2 /n.d.f. 

Df 

1.47(2) 

1.18(2) 

14.9(6) 

0.9 

Df 

1.66(3) 

1.77(4) 

36.2(36) 

1.0 

p) diag 

0.78(1) 

1.85(1) 

181.6(54) 

1.7 


TABLE II. Results of a fit of the maximally Abelian gluon 
propagators to (23 1 . 


B. Gluon propagator 

The gluon propagator in the maximally Abelian gauge 
has already been studied both in SU(2) JTSJQjaj and SU(3) 
[M US] for pure Yang-Mills theory. Here we extend 
those studies to full QCD with 2 + 1 flavors of dynam¬ 
ical quarks. In Fig. [2] we show the propagators of the 
coarse ensemble for the different quark masses including 
the quenched dataset. In addition to the three MAG 


C. Quark propagator 

We obtain the quark propagator in the standard way 
by inverting the (Asqtad) Dirac operator for a point 
source on a gauge field background. This is performed in 
Landau gauge and in the maximally Abelian gauge with 
C7(1) 3 x U( 1) 8 Landau residual gauge fixing. Addition¬ 
ally, we split the maximally Abelian gauge gluon fields 
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FIG. 2. The gluon propagator of the coarse ensembles of Table [I] renormalized at fj, = 4 GeV. 



FIG. 3. The gluon form factors of the fine ensemble. The 
solid lines illustrate the fit to Eq. (231. 


in their diagonal and off-diagonal components and invert 
the Dirac operator on the two parts separately. Thus, 
we obtain four “kinds” of quark propagators from each 
ensemble. 


1. Quark mass dependence 


On the five dynamical coarse ensembles of Table [TJ we 
calculate the quark propagator in Landau gauge, in max¬ 
imally Abelian gauge (MAG), on a pure diagonal MAG 
background (MAG-dia#) and a pure off-diagonal MAG 
background (MAG- off). The mass parameter of the va¬ 
lence quark propagator has been set to the value of the 
corresponding light sea quark mass for these five ensem¬ 
bles. This will allow for a systematic extrapolation to 
the chiral limit. 

In Fig. [4] the quark mass function M(k 2 ) is shown 
for the four types of gluon backgrounds (Landau gauge, 
MAG, MAG-diag and MAG-off) from all coarse ensem¬ 
bles, including a linear extrapolation to the chiral limit. 
In analog, Fig. [5] shows the corresponding quark wave- 
function renormalization functions Z{k 2 ). 

When comparing the Landau gauge quark propaga¬ 
tor to the quark propagator in MAG, the first obser¬ 
vation is that the MAG data, with the same statistics, 
results in more gauge noise from the Monte Carlo inte¬ 
gration. Moreover, the running masses of the MAG data 
lie higher than the corresponding Landau gauge masses, 
which holds from the largest mass of mi = 65.8 MeV 
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FIG. 4. The quark mass function from the coarse ensembles of Tabic [T| The bare quark masses are set to the values of the 
dynamical light quark masses. Additionally, a linear extrapolation to the chiral limit is shown. 


down to the chiral limit. The dynamically generated 
infrared masses, on the other hand, appear compatible 
within the error bars. 

It is evident that M(k 2 ) from the MAG -diag gluon 
background nicely resembles the Landau gauge analog 
(despite being more noisy), whereas M{k 2 ) from MAG- 
off gluons is constant, lying roughly 40% higher than 
the corresponding bare quark mass. In the chiral limit 
it is compatible with zero. Similarly, the wave-function 
renormalization function from MAG- off gluons comes 
out close to its tree-level value, Z{k 2 ) ss 1, independent 
of the quark mass. 

2. Infrared behavior 

The coarse ensembles cannot provide a clear picture 
of the infrared behavior of the quark propagators. To 
improve thereon we adopt the fine MILC ensemble of 
Table [TJ In order to keep the gauge noise and the simula¬ 
tion costs at an acceptable level, we use a valence quark 
with the mass of the heavier dyamical strange quark, 
m s = 68.0 MeV, instead of setting it to the light quark 
mass of m s = 6.8 MeV. 


In Fig. [6]we compare M(k 2 ) and Z(k 2 ) in MAG, MAG- 
diag and MAG- off, respectively, directly to the Landau 
gauge counterparts. While the MAG and Landau gauge 
mass functions agree within the error bars over the whole 
momentum range, Z{k 2 ) appears to be stronger IR sup¬ 
pressed in MAG as compared to Landau gauge. 

Similarly, the MAG -diag mass function agrees over the 
whole momentum range within the error bars with its 
Landau gauge counterpart. The renormalization func¬ 
tion exhibits qualitatively the same IR behavior as in 
Landau gauge. In contrast, the MAG -off quark dressing 
functions hardly show any non-trivial dynamics. This is 
the main finding of this study: it shows that the Abelian 
parts of the gluon fields not only dominate the purely 
gluonic interactions, but also the infrared interactions of 
quarks. 

V. SUMMARY 

We have fixed dynamical SU(3) lattice gauge fields 
to the combined maximally Abelian gauge and 17(1)3 x 
C/(l) 8 Landau gauge. From the lattice link variables we 
have extracted the continuum gluon fields which we sub- 
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|k| [GeV] |k| [GeV] 

FIG. 5. The quark function renormalization function from the coarse ensembles of Table [T] The bare quark masses are set to 
the values of the dynamical light quark masses. A linear extrapolation to the chiral limit is shown. All data have been been 
renormalized at /r = 3 GeV. 


sequently separated into purely diagonal (Abelian) and 
off-diagonal components. 

We investigated the gluon propagator from diagonal 
and off-diagonal gluon fields. Dynamical quarks lead 
to an IR suppression compared to the quenched case. 
The suppression becomes stronger when decreasing the 
quark mass. The screening is more pronounced in the 
MAG propagators compared to the Landau gauge propa¬ 
gator. Our findings confirm the manifestation of infrared 
Abelian dominance in the gluon propagator as found in 
earlier studies on quenched lattices. 

Finally, for the first time the maximally Abelian gauge 
quark propagator has been analyzed on a background of 
purely diagonal gluons as well as on the remaining, off- 


diagonal gluon background. The hypothesis of Abelian 
dominance implies that the non-Abelian gluon field does 
not propagate at a long-distance scale and hence that 
only the Abelian component is relevant at a long-distance 
scale. In accordance therewith, we have demonstrated 
that the quark propagator from a non-Abelian gluon 
background hardly shows any effects while the quark 
propagator from an Abelian gluon background closely re¬ 
sembles its Landau gauge counterpart. 
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